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SUMMARY 

Tho  report  covers  the  period  1  Ju.lv  197  2  -  30  June 
1973  which  is  characterized  by  continuous  imorove- 
ments  in  event  detection  and  location  performance 
during  routine  operation  of  the  array.  Tho  re¬ 
search  efforts  were  aimed  at  improving  the  event 
detectability,  and  the  potential  exploitation  of 
wave  scattering  effects  for  improving  UORSAR's  event 
detection  and  classification  capabilities. 

Work  completed  and  in  progress  is  presented  in  Chapter  ' 
II.  The  first  section  deals  with  seasonal  and  diurnal 
noise  level  fluctuations  plus  changes  in  the  character 
of  the  noise.  Next,  different  types  of  arrav  beam¬ 
forming  and  optimum  signal  estimation  techniques 
are  discussed.  The  importance  of  so-called  intrinsic 
time  and  amplitude  anomalies  or  wave  scattering 
(Cherno'')  effects  in  array  data  processing  are 
also  demonstrated.  Finally,  work  on  seismic  verifi¬ 
cation  problems  is  discussed  bxiefly. 


INTRODUCTION 

The  report  summarizes  the  NTNF/NORSAR  research  and 
development  efforts  during  the  interval  1  July  19  72  to 
30  June  1973.  In  the  first  part  of  the  period  most 
attention  was  given  to  development  work  like  software 
modifications  and  -implementation  of  new  data  processing 
routines.  For  example,  a  supplementary  event  detection 
processor,  based  on  so-called  incoherent  beamforming 
(Ringdal  et  al,  19  73)  ,  v/as  implemented  in  the  on-line 
system  in  September  1972.  For  the  purpose  of  editing 
a  daily  bulletin  of  seismic  events,  the  required  input 
data  may  be  read  directly  from  the  detection  log  tape. 

Thus,  a  daily  list  of  seismic  events  is  now  available 
every  morning  and  covering  the  previous  24  hours,  even 
if  the  Event  Processor  (EP)  is  behind  its  time 
schedule.  Moreover,  a  status  report  of  all  data  channels 
is  generated  daily  and  is  available  to  users  of  NORSAR  data. 
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The  research  topics  inves tiga ted  or  in  progress  are  mainly 
aimed  at  system  improvements  and  evaluation  of  the  array's 
event  detection  and  location  performance.  Also,  some 
aspects  of  the  event  classification  problem  have  been 
considered.  Most  attention  has  been  given  to  improving 
NORSAR's  event  detection  capability,  and  the  work  here 
comprised  optimum  amplitude  weighting  of  subarray  beam 
signals,  predictive  decomposition  models  for  explaining 
and  predicting  intrinsic  phase  shift  and  amplitude 
variations  across  the  array,  event  detector  false  alarm 
rate  fluctuations  and  signal-noise  wavelet  classification 
The  detectability  of  the  so-called  incoherent  event  detec¬ 
tor,  part  of  the  NORSAR  on-line  system,  is  superior  to 
that  of  conventional  beamforming  in  seismic  regions 
characterized  by  complex  P-signals.  An  evaluation  study 
of  NORSAR  event  detection  and  location  capabilities,  based 
on  the  array's  routine  performance  in  the  interval  Apr- 
Oct  ly 72  has  been  completed.  One  interesting  result  here 
is  that  the  NOAA-NORSAR  magnitude'  discrepancy  is  a  non¬ 
linear  function  of  magnitude,  i.e.,  NORSAR  reports  rela¬ 
tively  too  large  ra^  values  for  small  events.  Moreover, 
a  bias  analysis  of  NORSAR  event  magnitude  estimates  has 
also  been  undertaken.  Several  kinds  of  Vespugram  analyses 
of  core  precursor  phases  indicate  that  these  waves  are 
not  explainable  in  terms  of  the  standard  velocity  models 
for  the  earth's  core,  while  a  realistic  alternative  is 
scattering  sources  in  the  lower  mantle. 

The  above  research  topics  and  relevant  results  arc  described 
in  the  next  chapter.  In  general  only  the  main  results  are 
presented  here,  as  further  details  are  available  in  NORSAR 
Technical  Reports  or  from  papers  published  in  professional 
journals . 


3 


2  .  RESEARCH  LFFORTS 

The  research  activities  in  the  reporting  period  1  July 
1972  -  30  June  1973  have  been  focused  on  event  detection 
and  classification  problems.  Presently,  NORSAR  reports 
in  average  20  events  per  day,  but  recent  research  results 
indicate  that  significant  improvements  of  the  array's 
event  detectability  is  still  possible.  An  important 
problem  here  is  the  development  of  objective  criteria 
for  discriminating  between  weak  P-signals  and  signal 
shaped  noise  wavelets. 

Cone . tioned  on  the  present  NORSAR  computer  configuration, 
the  most  pressing  event  detection  problems  are  considered 
solved.  Thus,  at  the  end  of  the  reporting  period  more 
research  effort  could  be  spent  on  seismic  event  classifica 
tion  problems.  For  example,  conventional  discriminants 
criteria  have  been  adapted  to  NORSAR  data,  but  also  so- 
called  signal  space  expansion  techniques  are  under  in¬ 
vestigation.  In  the  reporting  period  a  number  of  visiting 
scientists  have  been  doing  research  at  the  NORSAR  data 
center  and  part  of  this  work  is  included  in  this  report. 

NORSAR  Event  Detec tor_Threshold_Se ttinq_and_ the 
False  Alarm_Rate 

At  NORSAR  a  significant  trend  in  seasonal  noise  level 
variations  occurs,  and  the  same  holds  on  a  diurnal  basis 
as  demonstrated  in  Figs.  1  and  2.  For  example,  extreme 
cases  with  a  variation  in  noise  power  up  to  18  dB  in  the 
frequency  band  2. 0-3.0  Hz  within  a  f cw  hours  have  been 
observed  at  a  large  number  of  NORSAR  short  period  sensor 
sites.  This  simply  moans  that  the  array's  event  detection 
capability  is  lower  during  winter  than  summer  and 
also  lower  during  the  day  than  the  night.  In  the  latter 
case,  there  are  roughly  20%  more  events  detected  during 
night  time.  Relevant  data  on  the  phenomena  are  presented 
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Bocim  average  of  LTA  for  the  time  period  6  Jan  -  2  3  Nov  1972. 
(LTA=Long  Term  Average,  which  is  equivalent  to  linear  pov.’er 
measured  in  a  window  of  approximately  30  sec.)  The  sampling 
rate  is  20  s/day  and  the  frequency  filtering  of  the  data 
from  which  the  LTA  is  computed  is  1.2- 3. 2  He. 


3  4 

TIME  OF  WEEK 


Fig.  2  LTA  .in  relative  units  as  a  function  of  time  of  week,  where 
day  1  is  Sunday.  (LTA  as  defined  in  Fig.  1)  Average  is 
made  over  46  weeks,  and  a  trend-removal  is  applied  to  the 
LTA  time  series  by  subtraction  of  daily  averages.  The 
sampling  rate  is  20  s/day,  and  the  frequency  fiJLcring 
of  the  data  from  which  the  LTA  is  computed  is  i.2-3.2  He. 
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and  discussed  in  a  recent  report'  by  Bungum  and  Ringdal  (1973)  . 
Also,  a  similar  study  for  long  period  noise  is  in  progress. 

An  important  but  mostly  ignored  aspect  of  noise  level  fluc¬ 
tuations  is  that  the  statistical  properties  of  the  background 
noise  change  too.  The  same  effect  is  also  obtained  by  using 
filters  with  different  passbands .  Henceforth,  from  the 
theoretical  studies  of  Rice  (1944)  and  Cartwright  and  Longue t- 
Higgins  (1956)  we  would  expect  that  the  noise  wave  train 
maxima  would  fluctuate  bctw?een  Gaussian  and  Rayleigh  prob¬ 
ability  density  distributions.  In  other  words ,  the  so- 
called  false  alarm  rate  would  vary,  i.c.,  the  number  of 
times  pure  noise  wavelets  trigger  the  event  detector  for  a 
fixed  SNR  threshold  (see  Fig.  3) .  This  phenomenon  does  not 
necessarily  mean  that  relatively  more  noise  wavelets  arc 
reported  as  seismic  events  under  adverse  noise  situations 


Fig.  3  Number  of  detections  as  a  function  of  SNR  for  two  time 
periods*,  covering  one  hour  of  day  time  and  one  hour  of 
night  time,  respectively.  The  dashed  line  has  a  slope 
of  -19.0. 


as  such  a  decision  rests  with  the  analyst.  Instead,  under 
favorable  conditions  too  few  events  would  be  reported  as 
the  SNR  threshold  value  in  the  event  detector  would  be  too 
large.  The  above  problem  was  first  considered  by  Lacoss  (1972) 
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who  forwarded  an  approximately  linear  relationship  between 
the  noise  stability  parameter  and  the  number  of  false  alarms. 
The  stability  parameter  was  defined  as  the  square  of  the 
noise  level  average  relative  to  its  variance.  Steinert  et 
al  (1973)  have  continued  this  work,  and  the  importance  of 
the  problem  and  also  the  potential  gain  by  having  a  float¬ 
ing  detector  threshold  setting  are  demonstrated  in  Fig.  4. 

Due  to  a  limited  data  base  the  SNR  threshold  values  range 
from  8  to  10  dB  in  the  figure,  while  the  corresponding 
values  in  the  operational  system  are  between  10  and  12  dB . 
Such  an  algorithm  for  the  prcspecif ied  false  alarm  rate 
is  feasible  to  implement  in  the  NORSAR  on-line  system,  and 
the  expected  gain  expressed  in  equivalent  SNR  units  would 
probcibly  amount  to  around  0.5  dB .  In  addition,  this  routine 
would  permit  a  more  efficient  computer  capacity  utilisation. 


co!sr  sv.ru  ty 


Fig.  4  False  alarm  rate  versus  noise  s  Labi  1.1  Ly  for  different  eve  nt 
detector  threshold  values.  Three  different  noise  situations 
were  analysed,  each  corresponding  to  one  he  ur  of  NORIAK  can- 
line  processing.  For  further  variation  of  the  noise 
structure,  three  different  bandpass  filters  were  also  used. 


5Y£!)]:_D2!'2£i:2£_k222d_on_Enyeloj2e_or  Incoherent  Beam  forming 

The  P-signals  recorded  by  NORSAR  are  only  partially  co¬ 
herent  across  the  array.  This  means  that  the  expected  gain 
in  sicjnal-to-noise  ratio  (SNR)  v/hich  is  proportional  to 
the  square  root  of  number  of  sensors  used  is  not  obtained 
during  array  beamforming  operations.  The  corresponding 
signal  energy  loss  increases  with  increasing  frequency, 
and  may  severely  degrade  the  array's  detectability  of  very 
short  period  P-waves.  This  problem  may  be  partly  circum¬ 
vented  by  replacing  or  supplementing  the  array  beam  traces 
with  the  average  of  subarray  beam  traces.  The  relative 
advantages  of  using  the  so-called  incoherent  beams  are 
modest  signal  losses,  better  estimates  of  the  noj.se 
variance  and  good  areal  coverage.  The  noise  suppression 
is  small  as  compared  to  array  beamforming,  but  could  partly 
be  compensated  for  by  using  high  frequency  bandpass  filter¬ 
ing.  As  mentioned  previously,  a  supplementary  event  detec¬ 
tor  based  on  incoherent  beams  was  implemented  in  the  on¬ 
line  system  in  September  1972.  Results  from  the  first 
two  months  of  parallel  operation  of  the  so-called  coherent 
and  incoherent  event  detectors  are  presented  in  Fig.  5 
and  Table  1.  The  improvement  in  the  array's  event  detect¬ 
ability  amounts  to  around  15  per  cent.  For  further 
details  see  the  report  by  Ringdal  et  al,  1972. 

The  characteristic  feature  of  a  seismic  array  is  real-time 
processing  of  data  from  a  large  number  of  sensors  organized 
in  a  certain  pattern  on  the  surface  of  the  earth.  As  is 
well  known,  when  sensor  separation  increases,  the  signal 
similarity,  in  general,  decreases.  The  consequence  here 
is  that  when  processing  signals  from  a  continental  array 
or  the  global,  seismological  network,  the  signal  suppression 
is  approximately  equal  to  the  noise  suppression,  resulting 
in  a  processing  gain  close  to  zero.  One  possible  way  to 
circumvent  this  problem  might  be  to  replace  the  individual 
signal  trace  by  its  envelope  as  we  intuitively  should 
expect  tills  kind  of  signals  to  exhibit  a  large  degree  of 
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signal  similarity  independent  of  sensor  separation,  seis¬ 
mometer  type,  etc.  This  hypothesis  has  been  tested  on 
WWSSN  station  records  (two  earthquakes  and  one  explosion) 
and  NORSAR  subarrav  beams  from  many  different  events.  The 
WWSSN  beam  pattern  for  an  earthquake  in  Chile  is  shov/n  in 
Fig.  6.  Signal  envelope  similarity  has  been  calculated 
through  cross-correlation  and  coherency  analysis.  Typical 
cross-correlation  values  v/ere  around  0.75  units  between 
WWSSN  envelope  signals.  Similar  results  were  obtained 
by  joint  analysis  of  22  different  NORSAR  events  in  the 
distance  range  3-145  deg.  Moreover,  using  data  on  the 
P-wavo  amplitude  variation  in  the  telescismic  distance 
range  and  the  theory  for  incoherent  event  detectors 
(Ringdal  et  al,  1972),  reliable  estimates  on  multiarray 
processing  gains  arc  obtainable.  It  is  interesting  to 
note  that  the  above  method  may  supplement  previously 
proposed  schemes  based  on  joint  detectability  analysis 
of  data  from  several  arrays.  The  above  topic  is  dis¬ 
cussed  in  some  detail  in  a  recent  report  by  Iluscbye  el. 
al,  197  2. 


Fig.  G  Response  pattern  for  the  Greeley  nuclear  explosion  in^ 
Nevada  12/20/1960  based  on  envelope  traces  for  19  WWSSN 
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lie  a  ml  or  mint}  or  simple  do  lay- and- sum  process'  g  is  exten¬ 
sively  used  in  analysis  of  P-wavos  recorded  by  the  large 
aperture  arrays  NORSAR  (Norway)  and  LASA  (Montana)  .  When 
tlit1  under  lying  assumptions  of  well— equ.ival  i.zcd  noise  levels 
and  identical  signals  between  instruments  .is  correct,  the 
corresponding  gain  in  SNR  is  optimum.  In  practice,  these 
restrictive  signal  and  noise  models  are  not  valid,  thus 
degrading  the  final  signal  estimate  and  the  event 
detectability  of  seismic  P-wavos. 

In  investigating  this  problem,  w c  followed  flic  line  of 
development  presented  by  Cliristof  fersson  and  Janson  (1973) 
where  the  interest  is  i. caused  on  the  relation  between  sig¬ 
nals  at  different  instruments,  i.c,,  the  space  spanned  by 
the  recorded  signals.  In  a  recent  paper  Cliristof  fersson 
and  ilusebye  (1973)  introduced  no  re  genem  alixed  P- signal 
models  during  array  bcamform.lng  and  the  corresponding 
least  squares  signal  estimation  techniques  'were  described. 
The  usefulness  of  these  data  processing  schemes  was  also 
demonstrated  in  analysis  of  more  than  one  hundred  LASA 
and  NORSAR  recorded  signals  (sec  Tables  2  and  3  and  Fig.  7) . 
lor  LASA  the  average  gain  in  SNR  relative  to  that  of  con¬ 
ventional  beamforming  for  the  telesoismic  event  was  around 
3.7  dB.  Most  of  this  was  obtained  by  accounting  for  noise 
level  vat iations  between  subarrays,  as  signal  coherency  aero: 
the  array  is  good.  For  NORSAR  the  corresponding  SNR  gain  wa: 
approx.  2.5  dB,  mostly  obtained  by  accounting  for  the  more 
complicated  signal  structures  in  this  case.  For  local 
events,  characterized  by  partly  incoherent  array  signors, 
relative  SNR  gains  amounting  to  5-10  dB  were  usually  ob¬ 
tained.  A  great  advantage  with  the  signal  estimation 
techniques  used  is  that  both  positive  and  negative  signal 
weights  are  permitted  thus  partly  avoiding  destructive 
signal  interference  during  beamforming  of  complicated  or 
very  weak  signals. 
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Very  weak  Wo a Lorn  Russia  event  recorded  by  IJORSAU.  A1 tentatively  the 
signals  may  represent  a  side  lobe  detection  of  a  local  explosion  in 
the  Haiti  c  Sea  area.  The*  columns  to  the  left  give  suharray  and  array 
beam  code*:,  (the  latter  equivalent  to  JJ  AM’. ,  JI J  ASH  in  Table  2)  ,  plot 
scaling  factors  and  relative  time  shifts  in  clr.oc .  The  r.igbthnnd  column 
gives  Model  1  IB  rubai  ray  weights,  SNR  and  relative  gain  in  dB  for  the 
array  beams.  Tin  signal  portion  marked  with  a  lino  is  6.5  see. 
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lDlrL'ii2£i£_^I^Sl72  Travel  Tine  and_Amulitude_An0mal2.es 
across_NORSAR 

As  is  well  known,  P-wave  travel  times  and  amplitudes  as 
observed  across  an  array  like  NORSAK  deviate  significantly 
from  that  expected  from  ray  theory  and  standard 
earth  models.  Except  for  special  subarray 
travel  time  correction  files  used  for  minimizing  signal 
energy  losses  during  beamforming ,  effects  of  the  above  types 
have  been  mostly  .ignored  in  mray  data  processing.  In 
recent  months,  considerable  efforts  have  been  spent  on 
analyzing  the  above  phenomena,  i . e . ,  whether  there  is  a 
non-random  pattern  in  the  P-v;ave  time  and  amplitude 
anomalies,  the  potential  .improvements  in  the  army  event 
detectability  and  classification  performance  by  tailing 
such  effects  into  account,  and  finally  to  find  more 
realistic  earth  models  to  explain  the  anomalous  P-wave 
propagation  effects.  The  results  obtained  so  far  will 
be  briefly  presented  in  the  next  subsections  starting 
with  time  anomalies. 

Trayol_Time_Anoma  1  i  c  s  _  a  c  r  o  s  s  _  t  h  e  _  N  0  R  R  A  R  _  A  r  ]■  a  y 

Array  beamforming  is  a  two-step  process J  first  the  in¬ 
dividual  subarray  beams  are  formed,  and  finally  the  array 
beam.  In  the  first  case,  the  necessary  time  delays  to 
ensure  proper  line-up  of  the  sensor  signals  are  based  on 
least  squares  P-wnvc  front  solutions.  In  array  beam- 
forming  the  P-wave  front  solution  used  is  a  first  order 
approximation,  while  the  second  order  terms  are  the  de¬ 
viations  between  observed  and  predicted  time  delays  (At^) 
using  so-called  master  events.  In  the  .latter  case, 
anomalies  amounting  to  tO.G  sec.  have  been  observed. 

This  justifies  an  analysis  of  the  effect  of 
ignoring  second  order  terms  in  subarray  beamforming. 
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The  results  obtained  give  that  the  subarray  travel  time 
anomalies  exhibit  a  distinct  regional  pattern  (see  Fig. 8) 
and  that  95  per  cent  of  the  At^  observations  have  values 
less  than  0.1  sec.  The  corresponding  subarray  bcamforming 
losses  are  around  0. 5-1.0  dB,  with  the  exception  of 
subarrays  0  5B  and  O'/C  (located  on  typical  Oslo  graben 
structures)  where  signal  energy  losses  may  amount  to 
2. 0-3.0  dB. 

Since  the  subarray  time  anomalies  are  non-random  quantities, 

they  should  lie  predictable,  so  the  following  experiment  was 

undertaken  (Dahle  et  al,  1973).  The  travel  time  variation tT .  , 

i 

across  NORSAR  is  modeled  as  a  function  of  two  factors,  namely, 
a  trend  effect,  equivalent  to  the  plane  wavefront  solution 

and  a  signal  or  wave  scattering  effect..  The  basic  idea  here 
is  physically  shown  in  Fig.  9,  and  the  corresponding  mathe¬ 
matical  formulation  is  given  in  eq .  (1). 

T.  =  T  +U  r.  +U  r.  TS.-tn.  (1) 

i  o  x  ix  y  ry  i  i 

where  r.  ,r.  arc  position  coordinates,  U,,U  are  trend 
ix'  ly  1  x  y 

components  (slowness),  S..  is  scattering  or  Chernov  (I960) 

effect,  and  n.  is  the  noise  at  the  i-th  site.  The  critical  factor 
l 

in  this  kind  of  analysis  is  the  autocovariance  functions 
typical  for  random  media  wave  scattering  models  (Chernov 
I960),  and  alternative1’/  observed  functional  values. 

To  demonstrate  the  usefulness  of  the  above  approach,  the 
arrival  times  at  roughly  half  of  the  NORSAR  sensors 
were  predicted  using  observed  travel  times  at  the  remain¬ 
ing  SP  instruments.  Using  actually  observed  travel  time 
data  as  a  reference  base,  most  of  the  intrinsic  travel 
time  anomalies  within  a  subarray  are  accounted  for  by 
inclusion  of  the  signal  effect  term  in  eq .  (1)  as 

demonstrated  in  Fig.  10.  Important,  the  minimum  of  the 
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Decomposition  model  of  coinprcnsional  scinmi  c  wave  field 
account/!  tuj  for  the  intrinsic  amplitude  and  travel  time 
fluctuations  observed  across  NORSAR  ro.i  smoiiK'tcr  sites. 
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sum  of  squared  differences  between  observed  and  predicted 
travel  times  was  observed  for  a  random  (Chernov)  medium  with 
correlation  distance  of  approx.  7  km  which  is  in  fair 
agreement  with  similar  results  obtained  for  LAS A  (Aki  197  >, 
Capon  197  3)  . 


P-wave  Ampl itude_Varia tion_across_NORSAR 

The  first  step  in  analysis  of  the  signal  amplitude 
variations  across  NORSAR  was  the.  probability  density  dis¬ 
tribution  of  this  parameter.  It  was  found  to  bo  approxi¬ 
mately  lognormal  (see  Fig.  11)  when  dominant  signal 
frequency  was  larger  than,  say,  0.9  II s.  This  result 
was  explained  by  Ringdal  et  al,  1972,  in  terms  of  multi 
plicative  response  effects  of  multilayered  earth  structures 
and  is  in  quantitative  agreement  with  the  wave  scattering 
theories  of  Chernov  (1960  )  and  Tartars); .i  ( 19 6  1)  . 


Fia  11  Observed  single  sensor  Amplitude  distribution  for  n  Kamchatka 
earthquake  occurring  Jan  03  at-06.3G.44  GMT  (NORSAR  bulletin). 
The  amplitude  values  wore  measured  after  applying  a  1.0- 3.4  Hz 
bandpass  filter.  The  normal  arid  lognormal  distribution  func¬ 
tions  estimated  from  the  observed  sample  mean  and  variance 
are  also  shown . 


- . »  ■  1 
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Seemingly,  the  P-wave  amplitude  variation  across  NORSAR 
is  random,  but  more  detailed  analysis  of  this  problem 
reveals  a  distinct  regional  dependent  pattern  in  the  in¬ 
dividual  sensor  or  subarray  beam  amplitudes  as  demonstrated 
in  Table  4. 


tablr  4 

Subarray  ranking  scores  for  different  seismic  regions.  In  all 
cases  the  results  obtained  are  significant.  For  details  on  the 
non-parametric  rank,  test,  see  Siegel  (1956). 
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The  prediction  experiment  of  NORSAR  observed  travel 
time  anomalies  mentioned  previously  was  repeated  on 
amplitude  data  and  some  results  are  shown  in  Fig  12. 

It  should  be  mentioned  that  the  Chernov 
media  parameters  used  in  modeling  the  autocovariance 
functions  and  giving  the  best  fit  to  the  observational 
data  are  similar  to  the  corresponding  values  ob¬ 
tained  in  the  time  anomaly  experiment. 

The  optimal  boamforming  procedure  (Chri stof fersson 
and  Ilusebye,  1D7  3)  discussed  in  a  previous  section 
actually  takes  advantage  of  the  skewness  in  the  sub¬ 
array  amplitude  distribution.  However,  this  method 
is  too  complex  for  on-line  data  processing,  but  a  viable 
alternative  is  using  the  simplified  scheme  of  masking 
the  weakest  subarrays  as  demonstrated  in  Fig.  13. 

Presently,  work  is  in  progress  to  map  the  NORSAR  sub  array 
amplitude  pattern  in  all  seismic  regions  expressly  for  tho 
purpose  of  'zero-one'  amplitude  weighting  during  on-line  array 
bcaruf orming .  The  subarray  cimplitudo  pattern  may  also 
be  instrumental  in  discriminating  between  very  weak 
seismic  signals  and  signal-shaped  noise  wavelets. 

Assuming  that  tho  above  optimal  boamforming  procedure 
is  used,  the  calculated  weights  should  be  matched  against 
that  expected  forget  given  region.  Preliminary  results 
give  that  this  method  would  be  an  important  diagnostic 
tool  in  classifying  weak  signals  -  noise  wavelets. 


DETECTABILITY  (dB)  DETECTABILITY  (<JB) 
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(WESTERN  RUSSIA) 


(GREECE- TURKEY) 


Fj.g.  13  Relative  gain  in  event  detectability  using  the  definition 
of  Kincjdal  et  al  (1972)  as  a  function  of  no.  of  NORSAR 
subarrays  for  different  bandpass  filters.  The  (a)  and  (c) 
figures  correspond  to  envelope  beamforming  using  (1)  1.6- 
3.6  Hz,  (2)  1.8-3. 8  Hz  and  (3)  2. 0-4.0  Hz  bandpass  filters. 
The  (b)  and  (d)  figures  correspond  to  conventional  beam- 
forming  using  (4)  1. 2-3.2  Hz,  (5)  1. 4-3.4  Hz,  and  (6)  1.6-3. 6 
bandpass  filters.  The  results  were  averaged  over  12  and  11 
events  respectively  for  the  'Western  Russia  and  Greecc-Turkoy 
regions.  All  the  events  analyzed  were  very  weak ,  i.e., 
having  SNR  values  between  2  and  4  and  thus  mostly  reported 
by  the  envfjlopo  detector. 


Hz 
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Seismic  Wave_Propaqa tion_ in_an_Earth_which_is_part ly 
Modeled  as_a_Random_or_Chernov_Media 

The  typical  features  of  the  Chernov  media 

are  small  perturbations  amounting  to  a  few  percent  of 

P  and  S  velocities,  density  and  the  elastic  parameters 

u  and  A.  The  corresponding  wave  scattering  effect 
1  x 

could  be  significant  as  shown  by  Haddon  (1973)  in  a 
theoretical  study  of  this  problem.  In  case  of  seismic 
arrays  this  kind  of  data  represents  an  excellent  tool 
for  observational  evidence  on  seismic  wave  scatter ing 
hypothesis  due  to  the  large  number  of  seismometers 
within  a  small  area.  For  example,  based  on  a  thorough 
analysis  of  NORSAR  recorded  core  precursor  waves, 
Doornbos  and  Ilusebye  (1972)  concluded  that  the  standard 
P-velocity  model  for  the  Earth's  core  probably  was  not 
quite  correct.  In  more  recent  studies  both  Doornbos 
and  Vlaar  (1973)  and  Haddon  (1973)  attributed  the 
above  precursor  waves  to  scattering  effects  in  the 
deep  mantle.  Moreover,  Aki  (1972)  and  Capon  (1973) 
have  used  array  data  for  mapping  the  extent  for  which 
the  crust  and  upper  mantle  beneath  LASA  could  be  con¬ 
sidered  a  Chernov  or  random  medium.  In  short,  based 
on  the  evidence  briefly  discussed  above,  we  feel  con¬ 
fident  that  wave  scattering  effects  could  be  used  as 
a  diagnostic  tool,  in  detecting  and  classifying  weak 
seismic  signals.  We  are  planning  to  investigate  most 
aspects  of  wave  scattering  effects,  partly  in  coopera¬ 
tion  with  Dr.  A.  Christof fersson ,  Uppsala  University, 
and  Dr.  R.A.W.  Haddon,  Sydney  University. 


Seismic_Magnitude_  Investigations' 

The  m,  magnitude  parameter,  measured  on  records  of  short 
b 

period  P-waves,  is  a  convenient  and  widely  used  tool 
for  ranking  of  earthquakes.  More  recently  this  parameter 


. . 


has  become  of  critical  importance  in  evaluating  event 
detection  and  discrimination  capabilities  of  various 
kinds  of  seismological  stations  and  networks.  The 
problem  of  a  possible  bias  in  the  NORSAR  estimation 
procedure  of  m^  magnitudes  and  also  that  used  by  the 
International  Seismological  Centre  (ISC)  in  Edinburgh 
have  been  investigated  by  Husebye  et  al,  1973.  The 
main  results  obtained  are  as  follows  (see  also  Table  5 
and  Fig .  14  )  . 


TABLE  5 


No .  of 

Subarrays 

dm ( loss) 

(m^-units) 

dm  ( skew) 

(m^-units) 

3 

0.28  1  0.06 

0.0  .<•  0.01 

6 

0.23  i  0.05 

-0.01  *  0.01 

9 

0.20  A  0.04 

i—i 

O 

o 

eg 

o 

o 

1 

12 

0.16  J  0.04 

-0.03  ±  0.02 

15 

0.14  1  0.04 

-0.04  '  0.02 

18 

0.11  t  0.03 

-0.05  +  0.03 

Operational 

0.08  I  0.03 

-0.07  i  0.03 

19<No<22 

_  _ _ 

Estimated  skewnef 
max.  power  distri 

s  of  subarray 
but  ion 

1.26  i  0.63 

Estimated  skewness  for  log- 
transformation  of  max.  power 

.-0.10  i  0.42 

Correlation  between  signal  loss 

-0.40  corr.  units 

and  skewness  effects 

Sample  size 

222  events 

Estimated  magnitude  biases  due  to  subarray  power  ]oss  and  skewed 
maximum  power  distribution,  conditioned  on  the  number  of  subarrays 


maximum  powe-  —  . 

The  latter  parameter  represents  a  decreasing  ordering  o  ^uinnay. 

based  on  naximum  power  ranking. 
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PREDICTED  MAGNITUDE  Irel.  scaie) 


Fig. 14  A  comparison  between  ovont  magnitudes  as  predicted  from  a 
multivariate  analysis  of  ISC  data  for  Japan,  and  that  of 
the  individual  stations  used  in  the  analysis.  The  relation¬ 
ship  betv/e  ' '  JSC  reported  magnitudes  and  predicted  magnitude 
is  also  given.  Dots  arc  observed  points  for  this  line.  This 
figure  .is  based  on  40  events  occurring  in  the  Japan  region 
in  1960  and  reported  jointly  by  the;  9  stations  listed  on 
the  figure. 

The  signal  energy  losses  observed  during  NORSAR  P-wave 
bearnf orming  do  not1  in  average  affect  its  event  magnitude 
estimates  due  to  a  skew,  approximately  lognormal,  P- 
amplitude  distribution  across  the  array.  A  comparison 
between  NORSAR-NOA/t  magnitude  gave  that  the  difference 
is  largest  at  m^~4 . 7  and  then  tapers  off  towards  both 
small  and  large  event  magnitudes.  A  multivariate  analysis 
of  ISC  data  for  Japan  and  the  /Aleutian  Islands  gave  a 
consistent  and  linear  relationship  between  the  ISC  event 
magnitude  and  that  predicted  from  subsets  of  5-9  sta¬ 
tions  in  the  rrij  4 . 0~6 . 0  magnitude  range  investigated. 

In  this  respect  the  ISC  reported  magnitudes  arc  con¬ 
sidered  unbiased.  We  also  found  that  the  magnitude 
observations  may  be  approximated  by  a  normal  distribu¬ 
tion.  In  many  cases  the  magnitude  station 
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correction  term  was  not  a  constant  but  a  function  of 
event  magnitude.  This  phenomenon  is  quantitatively 
explained  as  the  combined  effect  of  the  seismic  spectra 
scaling  law  (Aki,  1967,  1972)  and  the  crust-upper  mantle 
transfer  function. 

^D2iYsis_of_  the_0£c  rati  onal_Capabil.Ltins_  for_De  tec  tion 

§2^_LJ2£3£i2G_2f_§£i§1i1i£_Ey§!}i'S_at_WORSAR 

The  evaluation  of  the  NORSAR  event  detection  and  event 
location  capabilities  seems  to  be  a  popular  topic  as 
a  number  of  scientists  recently  have  worked  on  this 
problem,  namely,  Shlien  and  Toksoz  (1973),  Ringdal  and 
Whitelaw  (1973)  and  Bungum  and  Husebye  (1973).  The 
differences  in  the  results  presented  by  the  various 
authors  are  mainly  due  to  data  bases  covering  different 
time  intervais.  However,  as  Bungum  and  Husebye  (1973) 
used  both  the.  most  extensive  and  recent  data,  i.e., 
after  improved  time  correction  files,  better  bandpass 
filters,  incoherent  beamforming,  etc.,  had  been  incor¬ 
porated  in  the  array's  on-line  system,  we  prefer  to 
give  a  brief  summary  of  their  evaluation  of  NORSAR' s 
event  detection  and  location  capabilities  (see  also 
Table  6  and  Fig.  15  and  16) . 

Based  on  one  year  of  data,  Apr  1972-Mar  1973,  the  routine 
event  detectability  of  the  NORSAR  array  in  Norway  has 
been  investigated  in  terms  of  50%  and  90%  cumulative 
detectability  thresholds  which  were  derived  from 
frequency-magnitude  distributions.  The  best  performance 
was  observed  for  events  in  Central  Asia  and  adjacent 
regions  where  the  90%  cumulative  detectability  values 
are  in  the  range  3. 6-3. 8  NORSAR  m^  values.  For  tele- 
seismic  events  the  value  is  3..  8.  For  events  with  m. 
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Fig.  16  Average  difference  between  NOAA  and  NOR5AR  body  wave 
magnitudes  as  a  function  of  NOAA  maqnitude  for  all 
events  with  epicentral  range  30°-90°  and  110°-180° 
from  NORSAR  (regions  14  and  15  in  Table  6) .  The 
averaging  is  done  over  bands  of  0.3  magnitude  units, 
the  number  at  each  data  point  gives  the  number  of 
events,  and  the  upper  and  lower  bounds  are  the 
standard  deviations. 
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above  4.0  NOAA  reports  the  larger  value,  while 
NORSAR  reports  the  larger  for  events  below  4.0. 

The  accuracy  of  NORSAR-est imated  epicenter  solutions 
as  compared  to  those  of  NOAA  were  also  investigated. 

The  best  results  were  found  for  Japan  and  Central  Asia, 
where  the  median  location  difference  is  95  and  105  km, 
respectively.  For  telese.ismic  events,  the  value  is 
145  km.  The  biased  errors  in  the  location  estimates 
are  demonstrated  to  have  been  eliminated  for  most  of 
the  regions  considered.  Finally,  improvements  of 
the  present  NORSAR  event  detectability  performance  are 
discussed  in  view  of  recently  developed  array  data 
processing  techniques. 


Scismic_Verificat ion  Research 

So  far  our  main  research  efforts  have  been  aimed  at  im¬ 
proving  the  event  detectability  and  location  capability 
of  the  NORSAR  array.  Some  investigations  on  the  array's 
capability  to  discriminate  between  earthquakes  and  ex¬ 
plosions  have  already  been  undertaken,  although  at  the 
present  stage  only  conventional  class.if ^cation  criteria 
have  been  used  in  the  analysis  of  relevant  NORSAR  data. 


The  main  problem  with  the  application  of  ituj  :  criterion 

is  to  be  able  to  defect  the  surface  waves  from  smal 1  ex¬ 
plosions  and  earthquakes.  Three  major  signal  enhancement 
techniques  have  been  used  in  analysis  of  NORSAR  surface 
wave  data  with  good  results,  namely: 

-  Bandpass  filtering  centered  at  around  20  seconds, 
reducing  the  G  second  mi croseismic  energy,  which 
sometimes  can  be  very  s-trong  in  tiro  win!  or. 
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Beamforming,  which  works  well  if  the  noise  is 
separated  in  azimuth  from  the  signals.  The  signal 
similarity  is  always  high. 

Matched  filtering,  which  is  a  master  event  technique 
that  takes  advantage  of  the  time  invariance  of 
recorded  Rayleigh  waves. 

Fig.  17  shows  the  results  from  an  rn^  :  M  study  where 
those  techniques  have  been  applied  for  signal  enhancement. 
The  lowest  M  reported  is  2.5J  however,  at  other  times 
Mo  3.5  may  not  be  detected  due  to  variations  in  the 
background  noise.  Preliminary  results  from  an  m^  :  M , 
study  at  N0R3AR  have  been  published  by  Filson  and  Bungum 
(1972),  and  this  work  is  continuing. 


Fig.  17  Body  wave  magnitude  inj.,  versus  surface  wave  magnitude  M,. 

for  events  located  by  NORSAR  tn  Central  Asia  during  1971 
and  1972.  Open  c !r  Jes  indicate  presumed  nuclear  explosions. 


The  ni^  :  Mg  discrimination  criterion  v/orks  well  for 
larger  seismic  events.  However,  the  difficulty  of 
detecting  surface  waves  in  the  low  magnitude  range 
necessitates  investigations  of  event  classification 
capabilities  based  solely  on  P-waves .  Modified  versions 
of  the  complexity  and  third  moment  of  frequency 
criteria  have  been  tested  on  NORSAR  recorded  events  - 
earthquakes  located  in  Eurasia  and  North  America.  In 
the  latter  region  event  discrimination  using  P-waves 
only  is  relatively  poor,  and  also  the  array's  event 
detection  capability  is  not  specially  good.  On  the 
other  hand,  preliminary  results  indicate  that  fairly 
good  discrimination  between  earthquakes  and  presumed 
underground  explosions  is  achievable  for  Eurasia.  This 
statement  is  restricted  to  event  distances  larger  than 
around  3,000  km  from  NOUSAR,  as  the  modified  complexity 
criterion  does  not  give  satis! aetory  results  for  shorter 
distances.  The  principal  investigators  here  are  I. 
Noponen  and  1).  Rieber-Mohn. 


In  addition  to  the  seismic  noise  there  is  for  long 
period  waves  an  important  limiting  factor  for  the 
detectability  in  the  fact  that  waves  from  two  events 
are  very  often  interfering  with  each  other,  maybe 
as  much  as  20  per  cent  of  the  time.  The  long  period 
coda  from  a  large  event  may  last  for  hours,  and  another 
complicating  factor  is  that  the  energy  is  often  scattered 
in  azimuth  through  reflections  and  refractions  at  con¬ 
tinental  margins.  A  study  is  now  in  progress,  under¬ 
taken  by  H.  Bungum  and  J.  Capon  (M.I.T.  Lincoln  Lab), 
where  the  energy  distribution  in  the  coda  for  a  number 
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of  carefully  selected  events  is  studied  at  20  and  40 
second  periods.  The  advantage  of  working  at  40  second 
periods  is  that  the  multi  pathing  there  is  much  less 
severe  and  that  the  coda  fall  off  more  rapidly.  On 
the  other  hand,  some  events  may  have  energy  only  around 
20  second  periods.  The  results  for  NORSAR  are  comparable 
to  those  previously  obtained  for  I, ASA  by  Capon  (1.972), 
although  it  seems  that  NORSAR  data  gives  less  variation 
in  the  way  the  coda  around  40  second  wave  periods  fall 
off  with  time. 
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3 •  MISCELLANEOUS 

During  the  reporting  period  a  number  of  scientists, 
whose  names  are  listed  below,  have  visited  NORSAR  Data 
Processing  Center,  Kjcller,  for  special  research  purposes. 


D.  Doornbos 
Utrecht  University 
The  Netherlands 

I .  Noponen 

Se.ismological  Institute 
Helsinki,  Finland 

R. M.  Sheppard 
M.I.T.  Lincoln  Lab 
Cambridge,  Mass.,  U.S.A. 

II.  Ohlendorf 
Institut  fur  Geophysik 
Kiel,  West  Germany 

H.  Korhonen 
Oulu  University 
Fi nland 

S .  Pirhonen 

Seismological  Institute 
Helsinki,  Finland 

Professor  Tsujiura,  Tokyo, 

M.L.  MaKi 

Seismological  Institute 
Helsinki,  Finland 

A.  Christof fersson 
Statistical  Institute 
Uppsala,  Sweden 

E.  Iljortenberg 
Geodetic  Institute 
Copenhagen,  Denmark 

J .  Capon 

M.I.T.  Lincoln  Lab 
Cambridge,  Mass.,  U.S.A. 

J.  Vermeulen 
Utrecht  University 
The  Netherlands 


1  Jul  -  11  Sept  1972 
8  Dec  -  22  Dec  1972 
12  Jun  -  27  Aug  1973 

1  Jul  -  19  Dec  1972 
16  Feb  -  25  Feb  1973 
10  Jun  -  20  Jun  1973 

18  Sep  -  27  Oct  1972 


8  Sep  1972 


16  Nov  -  15  Dec  1972 


6  Nov  -  20  Dec  1972 

Japan  20  Dec  1972 

12  Dec  -  15  Dec  1972 

12  Feb  -  23  Feb  1973 

13  Jun  -  7  Jul  .1973 


21  Feb  -  14  Mar  1973 


16  May  -  4  Jun  1973 


18  Jun  -■  14  Sep  1973 
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NORSAR  scientists  participated  in  the  following  seminars, 
congresses  and  meetings  in  the  period  1  July  1972  -- 
30  Juno  1973. 


13th  General  Assembly  of  the  European  Seismological  Com¬ 
mission  in  Brasov,  Romania,  30  August  -  5  September  ,19  72. 
Participants:  K .  A .  Bertcusseri,  H.  Bungum,  II.  Gj0ystdal, 

and  E.S.  Ilusebye.  Altogether  the  NTNF/NORSAR  group  gave 
seven  talks,  wliich  are  listed  below: 

K.A.  Berteussen  and  E.S.  Ilusebye,  Seismicity  in 
terms  of  event  detection  thresholds 

II.  Bungum,  Event  detection  and  location  capabilities 
at  NORSAR 

H.  Bungum,  Array  stations  as  a  tool  for  microseismic 
research 

II.  Gj0ystdal,  E.S.  Ilusebye  and  D.  Rieber-Mohn, 
One-array  and  two-array  location  capabilities 
11.  Gj0ystdal  and  E.S.  Ilusebye,  Noise  suppression 
problems 

E. S.  Ilusebye  and  F.  Ringdal,  Multiarray  processing 
problems 

F.  Ringdal  and  E.S.  Ilusebye,  Event  detection  problems 
using  a  partially  coherent  array. 


Norwegian  Geophysical  Society  in  Ncsbyen,  Norway,  2-5  October. 
Participants:  11.  Bungum  and  E.S.  Ilusebye.  One  talk  was 

given . 


American  Geophysical  Union,  54th  Annual  meeting,  Washington, 
D.C.,  April  1973.  Participant:  K.A.  Berteussen.  One  talk, 
"Bias  analysis  of  NORSAR  and  ISC  reported  P-wave  magnitudes", 
by  K.A.  Berteussen,  A.  Dahle  and  E.S.  Ilusebye  was  presented. 


l'ouj.  Lh  Nordic  Seminar  on  Detection  Sc  is  mole  jy ,  Helsinki, 
Finland,  12-14  June.  Participants:  II.  Bungum,  a.  Dnhlc, 
H.  Gjtfystdal,  E.S.  Musebye,  N.  Maras,  D.  Riebor-Mohn , 

O.  Steinert.  The  NTNF/NORSAR  group  gave  ten  talks,  which 
are  listed  below: 


E.S.  Husebye,  A.  Dahlc  and  K.A.  Berteussen,  Analysis 
of  possible  non- random  errors  in  NORSAR  event  magnitude 
estimates 

D.  Riebcr-Mohn  and  I.  Noponcn,  New  short  period 
discrimination  criteria  used  on  NORSAR  events 

H .  Bungum  and  F.  Ringdul,  Diurnal  variation  of 
seismic  noise  and  its  effect  or  detectability 
0.  Steinert,  E.S.  Husebye  and  II.  Gj0ystdal,  Noise 
stability  and  false  alarm  rate  at  NORSAR 

E. S.  Husebye,  F.  Ringdal  and  J.  Fyen,  On-line  event 
detection  using  a  global  seismological  network 

II.  Gj0ys tdal ,  Array  detection  and  location  capabilities 
for  events  in  Central  Asia 

A.  Dahle ,  P-signal  variations  within  NORSAR  subarrays 

F.  Ringdal,  E.S.  Husebye  and  A.  Dahle,  Event  detection 
problems  using  a  partially  coherent  seismic  array 

A.  Chris toffersson  and  E.S.  Husebye,  Amplitude  weighting 
for  optimal  gain  in  SNR  during  array  beamforming 
0.  Steinert,  Stability  of  array  performance 


Norwegian  Geotravors  Meeting,  Bergen,  Norway,  4-5  May 
1973.  Participants:  K.A.  Berteussen,  II.  Bungum,  A.  Dahle, 
II .  Gj0ystdal,  E.S.  Husebye.  Three  talks  were  given. 
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